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Abstract 

The differential branching fraction with respect to the dimuon invariant mass squared, 
and the CP asymmetry of the —)• decay are measured for the first time. 

The CKM matrix elements \Vtd\ and |Us|; and the ratio |Ud/Us| are determined. 
The analysis is performed using proton-proton collision data corresponding to an 
integrated luminosity of 3.0 fb“^, collected by the LHCb experiment at centre-of- 
mass energies of 7 and 8 TeV. The total branching fraction and CP asymmetry of 
—)■ 7r^/r+/i decays are measured to be 

B{B^^ = (1.83 ± 0.24 ± 0.05) x 10“® and 

Acp{B^-^ = - 0.11 ± 0.12 ± 0 . 01 , 

where the first uncertainties are statistical and the second are systematic. These are 
the most precise measurements of these observables to date, and they are compatible 
with the predictions of the Standard Model. 
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1 Introduction 


The decay —)■ is a 6 —)■ d flavour-changing neutral-current process, which is 

suppressed in the Standard Model (SM)j^The suppression arises since the b —>■ 
transition proceeds only through amplitudes involving the electroweak loop (penguin and 
box) diagrams shown in Fig. 




Figure 1: Feynman diagrams of the penguin and box loop contributions to the b —)• di~^i process. 


In the SM, the top quark contribution dominates the loops, and an additional suppression 
occurs through the factor Vtd from the Cabbibo-Kobayashi-Maskawa (CKM) matrix. The 
decay is therefore sensitive to the presence of new particles that are predicted to exist in 
extensions of the SM, particularly in models where the flavour structure differs from that 
of the SM 00 The ratio of CKM matrix elements |ldd/kts| has been measured via 
and B^ mixing processes and b —>■ s((i )7 decays [^; it can also be determined 

from a measurement of the ratio of the branching fractions of the B~^ —)• decay to 

the more precisely measured B^ —decay 12 . Such ratios are also sensitive to 
the flavour structure of physics beyond the SM. 

The CP asymmetry of B^ —is defined as the relative difference between the 
decay widths, F, of the two charge conjugate modes. 


Acp = 


_ T{B —)■ 71 ) — F(i?+ —)■ 7r+/i+/i ) 


T{B- —)■ 7r“/x+/i“) + F(S+ —)■ 7r+/x+/i“) ’ 


( 1 ) 


The CP asymmetry is predicted to be non-zero due to interference between amplitudes 
that are proportional to the CKM matrix elements involved in the B~^ —)■ decay, 

namely VubV*^ and VtbV^^. Recent predictions for the CP asymmetry are given in Ref. [^. 
The B^ —)• decay was first observed by the LHCb collaboration 

branching fraction was measured to be 


13 and the total 


B{B~^ —)■ 7T~^^ ) = (2.3 ± 0.6 (stat) ±0.1 (syst)) x 10 ® . 


^Unless explicitly stated, the inclusion of charge-conjugate processes is implied. 
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This paper describes measurements of the differential branching fraction and CP asym¬ 
metry of the —)■ decay. The differential branching fraction is measured in bins of 

dilepton invariant mass squared, and normalised to B~^ —)■ J/'ip decays. These 

measurements are performed through fits to the invariant mass distributions. The branch¬ 
ing fraction and the ratio of the branching fractions B{B^ ^ 7r~^^~)/B{B~^ ^ 
are used to determine the CKM matrix elements \Vtd\ and and the ratio \VtdlVts\, 
respectively. The measurements are based on 3.0 fb~^ of pp collision data recorded using 
the LHCb detector at centre-of-mass energies of 7TeV and 8TeV. 


2 Detector and simulation 


The LHCb detector [14p5] is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing 6 or c quarks. The detector 
includes a high-precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4Tm, and three stations of silicon- 
strip detectors and straw drift tubes placed downstream of the magnet. The tracking 
system provides a measurement of the momentum, p, of charged particles with a relative 
uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The minimum 
distance of a track to a primary vertex, the impact parameter, is measured with a resolution 
of (15 -|- 29 /pt) where pt is the component of the momentum transverse to the beam, 
in GeV/c. The magnetic field polarity is inverted with a period of several weeks during 
data taking, which allows the charge asymmetries due to the detector geometry to be 
determined. 

The different types of charged hadrons are distinguished using information from two 
ring-imaging Cherenkov detectors. Photons, electrons and hadrons are identified by a 
calorimeter system consisting of scintillating-pad and preshower detectors, an electromag¬ 
netic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of 
alternating layers of iron and multiwire proportional chambers. The online event selection 
is performed by a trigger, which consists of a hardware stage, based on information from 
the calorimeter and muon systems, followed by a software stage, which reconstructs the 
full event. 

Samples of simulated 5+ —)■ 7r+/i+/i“, 5+ —and B~^ —)■ (/i+/i“)iL+ 

decays are produced from pp collisions generated using Pythia with a specific 
LHCb configuration 17 . Decays of hadronic particles are described by EvtGen [I^ , 
in which final-state radiation is generated using Photos (T^. The interaction of the 


generated particles with the detector, and its response, are implemented using the Geant4 
toolkit 20 as described in Ref. 21 . The simulated events are reweighted to account for 


known differences relative to the data in the transverse momentum spectrum of the B^ 
meson and the detector occupancy of the event. 
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3 Event selection 


Events are required to satisfy a hardware trigger, which selects muons with px > 1-48 GeV/c 
in the 7TeV data and px > 1.76GeV/c in the 8TeV data. In the subsequent software 
trigger, at least one of the hnal-state particles is required to have both px > 0.8GeV/c 
and impact parameter greater than 100 pm with respect to all primary pp interaction 
vertices (PVs) in the event. Finally, the tracks of at least two of the hnal-state particles are 
required to form a vertex that is signihcantly displaced from the PVs, and a multivariate 
algorithm is used to identify secondary vertices that are consistent with the decay of a 6 
hadron 
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Candidates are formed from pairs of well-reconstructed oppositely-charged tracks 
identihed as muons, combined with an additional track that is identihed as either a charged 
pion or a charged kaon for —)■ or —)■ decays, respectively. Each 

track is required to have a good ht quality, a low probability of overlapping with any other 
track, Px > 300MeV/c and to be inconsistent with originating from any PV. Candidates 
are required to have a good quality vertex ht and to be consistent with originating from a 
PV with the candidate’s momentum vector aligned with the direction between the primary 
and secondary vertices. 

Separation of the signal decay from combinatorial background is achieved using a 
multivariate classiher. A boosted decision tree (BDT) ^2 23 is trained using supervised 
learning with ten-fold cross validation 24 to achieve an unbiased classiher response. The 


background sample used to train the BDT consists of data from the upper sideband 
of the 7r ~^invariant mass distribution in the region greater than 5500 MeV/c^; the 
B~^ —)■ signal sample is obtained from the simulation. As no particle identihcation 

information is used in the classiher, it can be applied to both the pion and kaon modes. 
The features of the data that are used to classify the 7r+/r+/i“ candidate as signal- or 
background-like are the properties of the pion and muon tracks, and properties of the 


candidate. For the pion and muon tracks, the features used are the transverse 
momentum of the tracks, the impact parameter of the track, and the track quality. For 
the fjT candidate, the features used are the angle between its momentum vector and 
the direction vector between the primary vertex and the secondary vertex, and its hight 
distance, transverse momentum, and vertex quality. Two isolation variables and the 
absolute difference in momentum between each of the muons are also used in the classiher. 

The output of the multivariate classiher and the particle identihcation requirements 
are simultaneously optimised to maximise signal signihcance. Pseudo-datasets were 
constructed from simulated signal events and combinatorial background events taken from 
the upper mass sideband of data. Trial BDT and particle identihcation cuts were applied 
and an expected misidentihed-kaon component added to the pseudo-datasets. Wilks’ 
theorem 26 was used to determine a signal signihcance from hts to the pseudo-dataset. 


the value of which was passed to a maximisation algorithm that could vary the trial cut 
values. The classiher and particle identihcation cut values used to separate signal and 
background decays are chosen at the point of highest signihcance. Operating at this point, 
the classiher has a combinatorial background rejection of 99.8%, whilst retaining 66.9% of 
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signal events, and each event contains only a single candidate. As the classiher separates 
decays from combinatorial background, relatively pure samples of —)■ 

and )■ J/t/j fi~) K~^ events are also obtained using the same classiher reqnirements, 

when requiring a positively identihed kaon. 

The charmonium resonances are removed from the samples of B~^ —)■ 
and B^ ^ iji~ candidates by vetoing the regions 8.0 < < 11.0 GeVyc"^ and 

12.5 < < IS.OGeVyc'^. There are several other 6-hadron decays that could mimic the 

B~^ ^ signal. Decays such as 5^—)- Tr^vr^Tr''' and —)■ J/'il} , where 

there is double hadron-muon misidentihcation, are excluded from the B~^ —)■ 
dataset by muon identihcation criteria and the expected number of background events 
is found to be negligible. Partially reconstructed decays such as )■ , 

—)■ fi~ and —)■ fi~, where a kaon or a pion is missed, may 

satisfy the selection; however, simulation indicates that such events have a reconstructed 
mass that lies more than lOOMeV/c^ below the measured B^ mass. Therefore, such 
background events do not affect the signal yield extraction. 

There are two types of semileptonic decays that feature as backgrounds, 
B~^ ^ D^{K ^decays with kaon-muon misidentihcation, and the double semilep¬ 
tonic decay )■ D^{h^p,~V where can be a pion or kaon. The former decay 

is suppressed by requiring the to have a low probability of being a kaon. The latter 
decay has the same hnal state as the signal and cannot be completely removed by the 
selection. However, the distribution of double semileptonic decays as a function of the 
invariant mass varies smoothly, and can be modelled well in the ht from which 
the signal yield is extracted. The pion-kaon separation is not completely efficient: 6% of 
B^ —>■ p~ events are selected as B^ —)■ p~ events, and are modelled as a specihc 

background. The normalisation sample of B^ ^ J/'ijj {p'^p~)K~^ candidates is selected 
using the dilepton invariant-mass region around the J/ifj mass, i.e. 3096 ± 50MeV/c^. 
To remove much of the contribution from partially reconstructed decays, whilst keeping 
enough information to determine any effect on the signal, the 7r~^p~^p~ invariant-mass 
range 5040 < p^p~) < 6000MeV/c^ is used to extract the signal yield. 

4 Event yields 

The yields of 5+—)■ 7r+/i+/i“, 5+—)■ K~^p~^p~ and B~^ ^ J/'ip {p~^p~)K^ candidates are 
extracted by performing simultaneous, extended, unbinned maximum-likelihood hts to the 
invariant mass distributions m{'n'~^p'^p~) and m{K~^p'^p~) of the selected candidates. The 
total model for the invariant mass distribution is composed of a signal model, a combinato¬ 
rial background model, a model to describe partially reconstructed 6-meson decays and a 
model to describe 6-hadron decays with misidentified hnal-state particles. The signal model 
is an empirical function that consists of two Gaussian functions with power-law tails on 
both sides [^, and the same parameters are used for the B~^ ^ TT^p~^p~, )■ K~^p'^p~, 

and J/'ijj {p~^p~)K~^ decay modes. The model for the combinatorial background is 

described by a separate exponential function for each decay. In the B~^ —)■ p'^p~ data 
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sample, the misidentified —)■ decays where a kaon has been misidentihed as a 

pion, are described by a single Gaussian function with a power-law tail on the lower-mass 
side. The yield of misidentified —)• decays is constrained using the measured 

branching fraction 12 and the observed pion-kaon misidentification efficiency. The mass 
distribution of the misidentihed B~^ — fi^ii~ candidates is obtained by htting the 
invariant mass distribution of B^ —)■ J/ip candidates, where the kaon is required 

to have the pion mass, and which has been corrected to account for differences in the 
particle identihcation efficiencies that arise from the differing kinematics. The partially 
reconstructed B~^ decays in the 5+—)■ fi~ and the B^ ^ J/-!/) data are 

described by an empirical function, which consists of a rising exponential function that 
makes a smooth transition to a Gaussian function. This description allows the mixture of 
partially reconstructed 6-hadron decays to be limited to less than the maximum physical 
value of the B^ mass minus the pion mass, with a Gaussian resolution-smearing effect. 

The partially reconstructed 6-hadron decays in the B^ —sample are 
separated into three explicit components. Firstly, the double semileptonic decay 
B~^ ^ is included, as this is an irreducible background that ends at 

the B'^ mass. This is modelled by a falling exponential function that makes a smooth 
transition to a Gaussian function at high mass, where the parameters are hxed from a ht to 
simulated events. The yield of this component is left to vary in the ht. Secondly, the decays 
B~^ —)■ p+(7r+7r°)/i’''/i“ and B^ —)■ p°(7r+7r“)/i+/r“ are estimated to contribute a total of 
34±7 events to the data, from the measured branching fraction of B^ —)■ {7i~^7r ~)28 
and assuming isospin invariance. Lastly, the decay 5°—)- /o(7r+7r“)/i+/r“ is estimated to 
contribute 10 ± 2 events to the data, also below the B~^ mass. Each of these decays is 
modelled by a separate kernel-estimation probability density function (PDF) with a shape 
taken from simulated events reconstructed under the hypothesis. The yield of 

each of these decays has a Gaussian constraint applied with a central value and width set 
to the expected yield and its uncertainty. 

The invariant mass distributions of selected 7r+/i+/i“ and candidates are 

shown in Fig. along with the total htted model, signal component, and each background 
component. The ht gives yields of 94 ± 12 B^ —>■ 7T~^p~^p~, 2922 ± 55 B^ —)■ and 

(609.5 ± 0.8) X 10^ B~^ —)■ J/-^ candidates, where the uncertainties are statistical. 

The yield of B~^ —)■ 7r~^p'^p~ in each bin is given in Table [I] The ratio of GKM matrix 
elements is determined in the theoretically favourable bins 1.0 < g2 < 6.0GeV7c^ 
(low-g^) and 15.0 < < 22.0GeV7c^ (high-g^). The B^ ^ p~ yields are 879 ± 30 

in the low-g^ bin and 793 ± 28 in the high-g^ bin. The results of a simultaneous ht to the 
invariant mass distribution of B^ —)• p~ and B~ —)■ TT~p^p~ candidates are shown 
in Fig. and the measured yields are given in Table The small difference in total 
signal yield between this ht and that given in Table is due to the systematic effect of 
separating the background distributions by charge. Gonsistent results are obtained from 
datasets split between the two magnet polarities. 
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Figure 2: The fit to the invariant mass distribution of (left) selected )• candidates 

and (right) selected B+—)• candidates, with the total model and separate components 

as described in the legend. 


Table 1: The yields of ?■ decays in bins of dilepton invariant mass squared, with 

statistical uncertainties. 


bin (GeV7c^) 

B+^ 


' fl 

0.1 

-2.0 

22.5 

+ 

5.5 

4.8 

2.0 

-4.0 

7.5 

+ 

4.9 

4.0 

4.0 

-6.0 

11.1 

+ 

4.2 

3.5 

6.0 

-8.0 

9.5 

± 

3.9 

11.0 

- 12.5 

10.5 

± 

3.7 

15.0 

- 17.0 

9.7 

± 

3.3 

17.0 

- 19.0 

6.2 

± 

2.9 

19.0 

-22.0 

7.8 

± 

3.4 

22.0 

-25.0 

2.3 

+ 

2.1 

1.5 

0.0 

-25.0 

93.6 

± 

11.5 

1.0 

-6.0 

28.8 

+ 

6.7 

6.2 

15.0 

-22.0 

24.1 

+ 

6.0 

5.2 


The choice of models used for the partially reconstructed backgrounds, the semileptonic 
backgrounds, the misidentihed ii~ background, and the combinatorial background 

could all contribute as potential sources of systematic uncertainty. The dependence 
of the htted yields on these models is assessed by replacing the relevant component 
with an alternative model, as follows, and evaluating the change in yield in simulation 
studies and in the hts to data. The largest change in yield is assigned as the systematic 
uncertainty. Changing the models for the B’*'—)■ p+(7r’''7r°)/r ■*■/!“ and ^ 
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Table 2: The measured total yield from the simultaneous fit to the charge separated data, and 
the inferred yields of )• tt~^ fi~ and B~ —)• decays. 


J\f{B^ —)■ ) 

J\f{B+ -)■ Tr+n+fi~) 

Af{B —)■ TT /i+zi ) 

92.7 ± 11.5 

51.7 ±8.3 

41.1 ± 7.9 



Figure 3: The fit to the invariant mass distribution of (left) selected B~^ ^ 7r+/r+/i“ candidates 
and (right) selected B~ —)■ tt~ candidates, with the total model and separate components 
as described in the legend. 


decays to an exponential fnnction with a Ganssian high-mass endpoint contribntes 0.6% 
nncertainty to the measnred yield, and nsing an analogons shape for the 

)■ fi~ decays contribntes 0.7%. The parameters of the models are fixed to 

values obtained from a ht to the simulation. The systematic uncertainty of the model used 
for the semileptonic backgrounds is evaluated by allowing the exponent in the model to 
vary within the uncertainties produced by a ht to the simulation. This change contributes 
0.3% uncertainty to the measured B~^ —)■ 7r~^ yield. There is a negligible contribution 
from altering the model of the misidentihed decays or combinatorial background, and from 
changing the upper mass end-point of the ht range from 6000MeV/c^ to either 5500 or 
TOOOMeV/c^. 


5 Results 


5.1 Differential branching fraction 

The differential branching fraction of B^ in a bin of width is calculated 

relative to the normalisation channel >■ J/ip as 


di3(i?+—)■ 7r+/i+p ) 
dg2 


^B+-^iT+fi+fj,~ ■^B+—>‘J/'ilj{ii+ii~)K+ Ag^ 

( 2 ) 
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where M is the event yield, e is the total efficiency to select the decay, both of which 
are fnnctions of and )■ = (1.05 ± 0.05) x 10“^ is the measured 

branching fraction of the normalisation channel, with —?• = (5.961 ± 

0.033)% g. 

The total efficiency to select the candidates for the decays considered is computed from 
the product of the efficiencies to trigger, reconstruct and select the final-state particles and 
the candidate. This includes the geometrical acceptance of the LHCb detector and the 
efficiencies of the trigger and selection algorithms. These efficiencies are calculated using a 
combination of simulated signal events and data-driven methods. The use of the ratio of 
efficiencies of the decay modes ensures that many of the possible sources of systematic 
uncertainty largely cancel. The efficiency of the trigger depends on the kinematics of the 
muons, and this dependence contributes a source of systematic uncertainty relative to the 
signal yield at the level of 2%. The dependence of the particle identihcation efficiency 
on the kinematic distributions contributes a systematic uncertainty of < 0.1% for the 
muons, 2% for the pions and < 0.1% for the kaons. These uncertainties are evaluated 
by varying the binning of the kinematic variables, and include a contribution from the 
size of the calibration samples used. The calculation of the BDT efficiency is affected by 
small differences between the simulation and data. The dependence of the signal yield 
on these differences is assessed using the ^ J/-!/) (/r+/r“)ib+ and B~^ ^ J/'ip 
decays. The relatively large yield allows precise comparisons of data and simulation. 
The impact of using simulation to calculate the efficiency of the BDT is assessed using 
the observed differences between data and simulation in the normalisation channel; a 
systematic uncertainty of 1.4% is assigned. 

The measured values of the differential branching fraction are shown in Fig. |4 and 
given in Table The branching fraction agrees with SM predictions from Refs. 011 . 
although agreement in the lowest-g^ bin is only achieved when contributions from low-g^ 
resonances are taken into account, as in Ref. j^. The spectrum of candidates below 
1 GeV^c^ in a ±50 MeV window around the nominal B^ mass is shown in Fig. with 
hints of a peaking structure in the vicinity of the and u masses. The total branching 
fraction is computed from the integral over the measured bins multiplied by a scaling 
factor to account for the regions of not measured in this analysis. This factor is taken 
from simulation to be 1.333 ± 0.004, where the uncertainty combines the statistical and 
systematic uncertainties evaluated by using two different form factor models. The total 
branching fraction is therefore 

B{B^ —>■ Ti^jji^jjr) = (1.83 ± 0.24 (stat) ± 0.05 (syst)) x 10“® . 

The ratio of branching fractions of B{B^ ^ to B{B~^ ^ in the region 

1.0 < g2 < 6.0GeV7c^ is 

R(R+-)-7r+uV“) / N / N 

— ^ \ = 0.038 ± 0.009 stat ± 0.001 syst , 

o(i3+—)■ K + IJ,+IJi~) 
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•LHCb E3 APR13 M HKR15 ■ FNAL/MILC15 



(GeV2/c4) 


Figure 4: The differential branching fraction of >■ in bins of dilepton invariant mass 

squared, compared to SM predictions taken from Refs. 0 (APR13), (HKR15) and from 

lattice QCD calculations (FNAL/MILC15). 


and in the region 15.0 < < 22.0 GeV^c"^ is 

^\ = 0.037 ± 0.008 (stat) ± 0.001 (syst). 

0 ( 13 +—)■ 

These resnlts are the most precise measnrements of these qnantities to date. 

5.2 CKM matrix elements 

The ratio of CKM matrix elements |17d/f7s| can be calculated from the ratio of branching 
fractions, B{B^ —)■ —)■ 17+/r+/i“), and is given in terms of measured 

quantities 


\Vtd/Vts\^ = 


—)■ 7r+/i+/i ) / F^dg' 


X 


B{B+^ K+^+fi-) jF^dq 


(3) 


where is the combination of form factor, Wilson coefficients and phase space factor for 
the B~^ —)■ n{K) decay. The values of J F-^^xdq^ are calculated using the EOS package 29 , 
with B^ —)■ TT’*' form factors taken from Refs. 30,31 and 5+ —)■ 17+ form factors taken from 


Ref. 32 . The EOS package is a framework for calculating observables, with uncertainties, 
in semileptonic h-quark decays for both SM and new physics parameters. In order to 
take into account the correlations between the theory inputs for the matrix element ratio 
calculation, the EOS package is used to produce a PDE as a function of the B~^ —)■ 7r+/i+/i“ 
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Figure 5: The spectrum of ^ candidates in the region 0.1 - l.OGeV^/c^ in a 

±50 MeV window around the nominal mass, showing a peaking structure at 0.6GeV^/c^ 
that is in the region of the and u masses squared. 


and B~^ —)■ p, branching fractions in each of the relevant bins by Monte Carlo 

sampling of the theory nuisance parameters. A minimisation is performed to determine 

Table 3: The results for the differential branching fraction for B^ ^ p~ in bins of q^. The 
first uncertainties are statistical and the second are systematic. 


bin (GeV7c^) 

^{B+^n+p+p- 

) (10-9 GeV^c^ 

0.1 

- 2.0 

1.89 

+0.47 

-0.41 

±0.06 

2.0 

-4.0 

0.62 

+0.39 

-0.33 

±0.02 

4.0 

- 6.0 

0.85 

+0.32 

-0.27 

±0.02 

6.0 

- 8.0 

0.66 

+0.30 

-0.25 

±0.02 

11.0 

- 12.5 

0.88 

+0.34 

-0.29 

±0.03 

15.0 

- 17.0 

0.63 

+0.24 

-0.19 

±0.02 

17.0 

- 19.0 

0.41 

+0.21 

-0.17 

±0.01 

19.0 

- 22.0 

0.38 

+0.18 

-0.15 

±0.01 

22.0 

- 25.0 

0.14 

+0.13 

-0.09 

±0.01 

1.0 

- 6.0 

0.91 

+0.21 

-0.20 

±0.03 

15.0 

- 22.0 

0.47 

+0.12 

-0.10 

±0.01 
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\Vtd/Vts\, taking into account the data and this PDF, and the theory nuisance parameters 
are free to vary. The data are treated as uncorrelated between the two bins, but the 
full correlation between the theory parameters is accounted for. The value of the CKM 
matrix element ratio is determined to be 

— U.ZU_Q 04 , 


Vtd 

Vts 


where the uncertainty is the combination of the experimental (statistical and systematic), 
and theoretical uncertainties. Both contributions are approximately equal, and neither 
follows a Gaussian distribution. This is the most precise determination of |V)rf/V)s| in a 
decay that includes both penguin and box diagrams. 

Additionally, the values of \Vtd\ and |V)s| can be calculated via 


\Vtd 


2 




2 


—)■ TT+zi+zi ) 

B{B+^ ib+ZrV") 


and 


/F^dg2 


(4) 

(5) 


where EOS is used to compute the theoretical input. Combining the results from the high- 
and low-g^ bins gives 


\Vtd\ = 7.2;o| X 10-3 and 
|Gi,| = 3.2;°;^xlO-2, 


where the uncertainties are due to both the branching fraction measurements and the 
theory nuisance parameters. As the |Gd/hts| determination uses both the S’*"—)■ tt+z^+z^” 
and B^ -2- branching fraction measurements, the theory nuisance parameters 

take different values to those in the separate \Vtd\ and \Vts\ determinations, where only 
one of the branching fractions is used. The ratio of |V)d| and \Vts\ is therefore not identical 
to the measurement of given above. The uncertainty on \Vtd\ has approximately 

equal contributions from experimental and theoretical uncertainties, while the uncertainty 
on |V)s| is dominated by the theoretical uncertainty. 


5.3 CP asymmetry 


The CP asymmetry of B'^ -2- vr^z^+z^ , as dehned by Eq. can be computed from the raw 
yield asymmetry. 


_ N'{B -2- TT fi ) — J\f{B^ —)■ TT+Zi+Z^ ) 

RAW j\f(^B- —)■ 7r“Zi+zi-) +N'{B+ —)■ TT+zi+z^-) ’ 


( 6 ) 


where M is the signal yield for the given decay-mode. This raw asymmetry is corrected for 
the production asymmetry of the B^ mesons and the detection asymmetry of the decay 
products, under the approximation 


•AcpiB'^ —)■ ) — .Araw — .Ap — A.det , 


(7) 
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where is the i?^-meson production asymmetry, and .Adet is the detector asymmetry 
for the pions and muons. 

The production asymmetry of and B~ mesons at LHCb has been measnred to 
be (—0.6 ± 0.6)% nsing the B^ ^ J/'ip decay |33 . The momentu m spectrnm 

differences between the B'^ —)■ J/'i!) iJi~)K^ and B^ -x 7r+/i+/i“ decays are fonnd to 
have a negligible impact on this asymmetry. The charge asymmetry of the LHCb detector 
for 7r+ and tt~ has been measnred in decays to be = 0.9914 ± 0.0040 

and Et^+Iet^- = 1.0045 ± 0.0034 for the two magnet polarities. These efficiency ratios 
give detector asymmetries of (—0.43 ± 0.20)% and (0.22 ± 0.17)% for the two magnet 
polarities, where the differences in the momentnm spectrnm are acconnted for in bins of 
momentum, transverse momentum and azimuthal angle. The relative tracking efficiency of 
differently charged pions is consistent with unity when averaged over the the two magnet 
polarities 34 . The pion identihcation asymmetry is derived using —)■ K~ti^ decays and 


is calculated to be less than 0.087% when momentum spectrum differences are accounted 
for. Additional effects from the prodnction and detection asymmetries are negligible and 
do not contribnte to the hnal systematic nncertainty. 

The raw CP asymmetry, AIraw, of the B^ -x iJi~ candidates is measnred to be 
—0.11 ± 0.12. The value of Acp for B^ ^ is calculated to be 


AcpiB 




/i+/i ) = 


-0.11 ± 0.12 (stat) ± 0.01 (syst), 


which is consistent with a recent SM prediction 1^. 


6 Summary 


A measnrement of the differential branching fraction of the decay B~^ -x has 

been presented, and is found to be consistent with SM predictions, and to have a possible 
contribntion from B~^ —)■ p^{u)7i~^ decays. The CP asymmetry of the decay has been 


measnred and is consistent with a recent SM prediction [^|. The valnes for the CKM 
matrix elements |T)d| and |V)s|, and the ratio \Vtd/Vts\ have also been determined, and 
are in agreement with previons measnrements. These results constitute the most precise 
measurements to date of a 6 —)■ transition and supersede those of Ref. 13 . 
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